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ABSTRACT 

The Cold Dark Matter (CDM) model has an elegant simplicitly which makes it very 
predictive, but when its parameters are fixed at their 'canonical' values its predictions 
are in conflict with observational data. There is, however, much leeway in the initial 
conditions within the CDM framework. We advocate a re-examination of the CDM 
model, taking into account modest variation of parameters from their canonical values. 
We find that CDM models with n = 0.8-0.9 and h = 0.45-0.50 can fit the available 
data. Our "best fit" CDM model has n = 0.9, h = 0.45 and Cj /C^ = 0.7. We discuss 
the current state of observations which could definitely rule out this model. 
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NTRODUCTION 

"Qp most successful of all cosmological models for large-scale 
^tljicture has surely been that of inflationary (fitot = 1, ini- 
t^ly scale-invariant fluctuations) cold dark matter, other- 
wise known as CDM (Peebles 1982; Blumenthal et al. 1984; 
Davis et al. 1985). It has been thoroughly examined via N- 
body and SPH simulations on scales ranging from galaxies to 
clusters and superclusters. Only in the post- COBE (Smoot 
et al. 1992) era, when the normalization of the fluctuations 
in the linear regime has finally been determined, has the the- 
ory become vulnerable to attack on several fronts. The over- 
whelming opinion among cosmologists currently holds that 
CDM is now in flagrant conflict with observational data, 
notably on scales of ~ 1 — 10 Mpc (e.g. (Davis et al. 1992; 
Liddle & Lyth 1993; Ostriker 1993)). With the popular al- 
ternatives of either a cosmological constant-dominated or an 
open universe requiring initial conditions that are (for many 
physicists) distasteful to contemplate, we think it behooves 
us to carefully reconsider the CDM model before discarding 
it. 

Indeed, we will show that remaining within the frame- 
work of the standard CDM model but with judicious tuning 
of the canonical parameters, cold dark matter is worthy of 
serious reconsideration. New observations will be needed, 
described below, before a definitive resolution of the astro- 
nomical issues can be accomplished. We note in passing that 
confirmation of a small (~ few eV) neutrino mass would 



provide another means of maintaining a viable cold dark 
matter-dominated universe, in this case with a ~ 20% ad- 
mixture of hot dark matter. However in the absence of an 
experimental measurement of neutrino mass, we regard the 
coincidence required in order of magnitude between and 
f^CDM to be unappealing. In any case we feel that a mixed 
dark matter model is best considered as a minor variation 
on the theme of a CDM universe. 

We will consider models which keep the attractive fea- 
tures of standard CDM, while incorporating only mild devi- 
ations from the standard parameters. The main attraction 
of CDM is its simplicity: (1) Qtot = 1, with the dominant 
mass in the form of massive particles whose interactions and 
velocity dispersion are negligible; (2) adiabatic fluctuations 
drawn from a Gaussian random distribution; and (3) initial 
conditions which can arise naturally from inflation. We be- 
lieve that there are models which satisfy these requirements 
while being in reasonable agreement with the data. 



2 COSMOLOGICAL DATA 

Given the enormous breadth of available cosmological data, 
which pieces of information does the CDM model need to 
address? Making such a choice is inevitably a subjective 
process that involves taking a somewhat skeptical view of 
many published claims, in part because of uncertain sys- 
tematics such as those encountered in recent Hubble con- 
stant determinations (cf. (Freedman et al. 1994; Nugent et 
al. 1995)). We will also avoid the issue of whether there is 
any evidence for or against having Qo = 1 on the largest 
scales as measured by velocity flows (Dekel 1994; Strauss 
& Willick 1995; Shaya, Peebles & TuIIy 1995); if it be- 
comes clear that Qo < 1 then the models we discuss will 
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no longer be viable. Instead of these classical cosmologi- 
cal parameters, we will use data which measure the power 
over some range of scales: (T(r) is the rms overdensity within 
spheres of size r, and the dimensionless power spectrum is 
A^(A;) = dal^/dlnk = P {k) I {2t:^ ) . 

We will specifically concentrate on the following list of 
constraints: 

(i) the COBE normalization at the largest scales; 

(ii) a measurement of the shape of the spectrum at large 
scales; 

(iii) determination of the normalization of the fluctua- 
tions on cluster scales; 

(iv) the potential energy as measured by velocities on 
Mpc scales; and 

(v) an indication that there may be an acoustic peak in 
the CMB fluctuations. 

To commence, we remark that the theoretically cleanest 
number is the COBE (Bennet et al. 1994) normalization of 
the potential fluctuations on the horizon scale (Bunn, Scott 
& White 1995; Stompor, Gorski & Banday 1995). We impose 
the COBE constraint by normalizing all of our models di- 
rectly to the COBE 2-year data using the method of (White 
& Bunn 1995). While the COBE data gives information on 
the shape of the power spectrum at large scales, the change 
in likelihood across the range of parameters we consider is 
not significant, so we do not include it. 

Secondly we want to use information from large scales 
about the break in the spectrum. Frequently this is specified 
in terms of a shape parameter, F ~ Qo A, but we prefer to use 
a measure more closely tied to the data. Since information 
on scales > 100 /t"^ Mpc is sparse, most current F limits 
are really statements about relative amounts of "large" to 
small scale power. We choose to consider the ratio between 
the contribution to the rms power at large (~ 30 h'^Mpc) 
scales, Ao.o5 = A(A; = 0.05 AMpc""'), to that at small scales: 
ag (the rms overdensity at 8/t~^Mpc, which probes power 
at wavenumber ~ 0.2/t Mpc""' )* both in linear theory. A 
value for Ao.os /f 8 can be found from the data compilation of 
(Peacock & Dodds 1994): Aom/(Ts = 0.27±0.02. This value 
for Ao.os/fs is consistent with the value determined from 
power spectral analysis of a range of surveys once corrections 
are made for redshift space distortions. For standard CDM 
this ratio is 0.22, which is the usual statement that this 
model has insufficient large-scale power. If we fix ra = 1 and 
use the fitting function for T{k) of (Efstathiou 1990), we find 
Ao^o5/(T8 ~ 0.21 -I- 0.22(1/2 -F) for F near 1/2- While there is 
certainly more information about large-scale structure than 
this, in general a smooth P{k) that passes this test, and 
others outlined below, provides a reasonable fit to most of 
the data (see Figure 2). 

We next need to confront the overall amplitude of the 
mass fluctuations on cluster scales: ag. The best way to es- 
timate this is probably by means of cluster abundances, 
which can be related to the amplitude of the fluctuations 
in linear theory through the Press-Schechter ansatz. Recent 
determinations (Bond & Myers 1991; White, Efstathiou & 

* This may be compared with the "excess power" defined in 
(Wright et al. 1992), which probes similar scales but is harder 
to obtain from the published data. 



Frenk 1993; Carlberg et al. 1994) lead to as ~ 0.5 - 0.8, 
where the principle uncertainty is in the assumed cluster 
mass ((78 ~ m" * in the Press-Schechter theory), though the 
formalism employed and data used give rise to some of the 
spread. A lower limit on ag can be obtained by considering 
rare (high mass or temperature) clusters at redshift z > 0, 
though this is not a constraint for the models considered 
here. Since for optical galaxies ((78, g)^ ~ 0.9 ±0.05 (Foveday 
et al. 1992; Foveday et al. 1995), this constraint can also be 
phrased in terms of bias: b = ag^g/ag ~ 1.2-2. 

One area where the CDM model has run into trou- 
ble has been its prediction of velocities on Mpc scales 
which are too "hot" (Suto, Cen & Ostriker 1992; Gelb & 
Bertschinger 1994; Schlegel et al. 1994; Nolthenius, Klypin 
& Primack 1994). This is a problem for most Qo = 1 
models, which predict hot velocity fields whereas the small 
scale velocity field appears locally to be quite cold (Pee- 
bles 1992). Rather than focus on the relative pair dispersion, 
(712, which is dominated by clusters, we shall consider the 
small scale rms velocity of individual galaxies: ?)rms. Red- 
shift catalogues typically give 3-dimensional measurements 
of Vrrns ~ 350 km s"^ (Miller & Davis 1995), though refined 
analyses of N-body simulations and large redshift catalogues 
will be required for more precise constraints. To address this 
we calculate the 2nd moment of the mass correlation func- 
tion, i.e. J2(l /t"^Mpc), in non-linear theory using the fitting 
function for P{k) of (Mo, Jain & White 1995). This gives an 
estimate of the potentials in a filtered version of the Fayzer- 
Irvine or "cosmic energy" equation, which can be related to 
the three-dimensional rms velocity on Mpc scales through 
(Peebles 1980; Davis & Peebles 1983) 



Note that since conventionally J2 (1 A""' Mpc) is measured 
in (/t~^Mpc)^, the dependence on Hq cancels in Viras- On 
purely dimensional grounds J2(l/t~^Mpc) should also be 
related to the pairwise velocity dispersion squared, (7^2, de- 
termined by the cosmic virial theorem, if scaling solutions 
are used to relate the three-point function to the two-point 
function (as is often done). We note the accurate calcula- 
tion of the pairwise velocity dispersion is not unambiguous 
(e.g. (Zurek et al. 1994; Marzke et al. 1995)). Given the 
above, it seems that a reduction in J2 on Mpc scales by a 
factor > 3 (relative to standard CDM) is desirable. We find 
that this it is possible to reduce J2 by a factor 3-4, but it is 
unclear whether this is enough. 

At essentially the same scale, but constraining the 
power spectrum from the other direction, are abundances 
of quasars and damped Lya systems at high redshift. We 
can regard the data as leading to a lower limit on the rms 
mass fluctuation at 1 /t~^Mpc in linear theory: (7i. This can 
be used to compute abundances based on Press-Schechter 
(or peak-patch) estimates and to obtain a feel for the red- 
shift at which structures would form: (1 -|- Zform) ~ (7i (but 
see (Katz et al. 1994) for simulations of early structure for- 
mation). Ideally we would like to keep (7i > 3, though we 
shall see that early structure formation is not a problem for 
any of our models. 

Finally, there is the observationally immature but the- 
oretically crucial issue of the acoustic peak in degree-scale 
CMB fluctuations (Scott, Silk & White 1995). Although the 
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possibility of systematic errors and foreground contamina- 
tion remain an issue, the data at face value require the ratio 
of the large angle to degree-scale fluctuations -D200 /-Dio ^ 2 
(with a rather large error bar) where Di ~ i^Ci. This last 
constraint will be very important if strengthened, since the 
peak goes away if we tilt too much or introduce too much 
tensor contribution. However, this depends to a large extent 
on Q_b: the peak becomes higher as Q_b increases. If new 
measurements of the primordial deuterium abundance give 
a value of a few xfO~* (Songaila et al f994), forcing Q_b to 
the lower end of the nucleosynthesis range, then the presence 
of a Doppler peak would put a much stronger constraint on 
the allowed n and hence ag . 

Should a variant of CDM be able to fit these numbers, it 
will simultaneously fit many of the other less rigorous tests 
summarized elsewhere (e.g. (Silk 1994; Bond 1994; Bah- 
call 1994)). We regard the strongest of these additional tests 
as being the requirement that the Hubble constant should 
not be much lower than Ho ~ 50kms~^ Mpc~^, and that 
the age of the universe should be > 12 Gyr. It appears to be 
possible to match all of the above constraints by considering 
small deviations from the parameters of standard CDM. 



3 THE CDM MODEL 

In almost every 'realistic' inflationary model, the spectral 
index, n, of the primordial fluctuation spectrum is only 
approximately unity. Indeed (unless there is tuning of pa- 
rameters e.g. (Barrow & Liddle 1993)) the limit of exact 
scale invariance corresponds to an exactly de Sitter infla- 
tionary phase (which would never end) and is unphysical. 
This is apparent in the divergence of the amplitude of the 
adiabatic fluctuations generated by inflation as ra 1. 
Two classes of inflationary models generate fluctuation spec- 
tra which are almost exact power laws: natural inflation 
(Adams et al. 1993) and power law or extended inflation 
(Lucchin & Matarrese 1985; La & Steinhardt 1989; Stein- 
hardt & Accetta 1990; Kolb, Salopek & Turner 1990; Lyth 
& Stewart 1992). The other classes generally have spectra 
which are only approximately power laws, with the correc- 
tions becoming more important as the spectrum departs 
from scale invariance (see e.g. (Stewart & Lyth 1994)). 
For simplicity, we will focus on power law initial fluctua- 
tion spectra, though we emphasize that a scale-dependent 
spectral index is not unlikely if the fluctuations are gen- 
erated by inflation (Steinhardt & Turner 1984; Kofman & 
Linde 1987; Kofman & Pogosyan 1988; Salopek, Bond & 
Bardeen 1989; Liddle & Turner 1994). Many authors have 
stressed that having n somewhat below unity is a desir- 
able feature (e.g. (Bond 1994; Scott, Silk & White 1995)), 
and certainly tilted CDM models have been considered be- 
fore (e.g. (Cen et al. 1995; Liddle et al. 1992; Muciaccia et 
al. 1993)). However such studies have usually been in terms 
of relatively extreme tilt, i.e. n < 0.7, and without using 
the accurate normalization on large scales provided by the 
COBE 2-year data. 

Power law inflation models predict a component of ten- 
sor fluctuations whose amplitude is related to that of the 
scalar component by an amount proportional to (1 — ra) (Lid- 
dle & Lyth 1992; Davis et al. 1992; Stewart & Lyth 1994), 
while natural inflation models predict a negligibly small ten- 



sor component (Adams et al. 1993). The tensor fluctuations 
affect only the predictions of the model for large scale CMB 
anisotropics, but this affects the normalization of the model 
in terms of COBE. The ratio of tensor and scalar compo- 
nents is usually expressed in terms of their contribution to 
the quadrupole CMB anisotropy. For power law inflation 
Cf/Cf = 7(1 - ra) while for natural inflation Cf x 0. We 
examine tilted CDM models both with and without tensor 
contributions and restrict the tensor contribution to that 
predicted by power law inflation. 

We attempt to fit a wide range of observational data by 
allowing ra, Cj/C;f and h to deviate by only small amounts 
from their canonical values. Specifically, we assume that 
QbA^ = 0.02 from BBN (Copi et al. 1995; Krauss & Ker- 
nan 1995), and allow h to range between 0.45 and 0.55 and ra 
to range between 0.8 and 1. Higher values of ra and h could 
clearly be considered, but these give worse and worse fits 
to the data. Adopting a lower value of Q_b would make it 
more difficult to keep -D200/-D10 ^ 2, since we are generally 
decreasing ra to get ag to observationally acceptable levels. 
Though in principle it would also become harder to fit the 
high redshift damped Ly a abundance with lower values of 
fis, our models produce enough small scale power that this 
is not a serious constraint. 

Power spectra for our models were calculated numer- 
ically with a Boltzmann code. The transfer functions, in- 
cluding the effects of our assumed baryon content, can be 
fit by the form of (Efstathiou 1990) with F = 0.39, 0.45, 
0.51, for h = 0.45, 0.50, 0.55 respectively^. Fitting func- 
tions for various parameters discussed above, over the range 
0.45 <h< 0.55 and 0.8 < ra < 1, are given in Table 1. The 
error on ag and ai simply from the COBE normalization 
is 7.5%, while it is 15% for J2 (1 A""' Mpc). The expressions 
in Table 1 are formally accurate to < 5%, but for J2 have 
additional uncertainty due to the non-linear nature of the 
calculation. 



4 DISCUSSION 

Figure 1 shows the contours of Ao.os/fs, fs, fi and 
^2(1 /i~^Mpc) vs ra and h for models without and with ten- 
sors, respectively. To have ag < 0.8 we are forced into the 
bottom left region of the plots. Choosing round numbers 
our best fit models have either: (1) h ~ 0.45, ra ~ 0.8, 
C2^/C2^ ~ 0; or (2) h ~ 0.45, ra ~ 0.9, C2^/C2^ ~ 0.7. 
The values from Figure 1 are listed in Table 2 for these 
two models and for "standard" CDM. Also shown in Fig- 
ure 2 are the matter power spectra for these two mod- 
els plus an MDM model (with = 0.2) and a ACDM 
model (with Qq = 1 - f^A = 0.3, h = 0.8, ra = 0.95 and 
C2^/C2^ = 7(1 - ra) (Scott, Silk & White 1995)) for compari- 
son. The data points are an estimate of the linear P{k) from 
(Peacock & Dodds 1994), and are plotted assuming Qo = 1. 
[The non-linear power spectrum we have used to estimate 
J2 (1 /i""' Mpc) and v^ras is also shown.] There is an uncer- 
tainty of at least 20% in the overall normalization of these 
data points. In addition the amplitude of the points scales as 

^ Note: for ri < 1 the "effective" P for P{k) would be lower than 
for T(k). 
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0(7" '^, so we have multiplied the ACDM curve by 0.3" '^ to 
compensate. The radiation power spectra, Di = £{£ + l)Ci 
vs £, for our two models are shown in Figure 3. 

Our models are able to accomodate the measure of 
Ao.os/fs from (Peacock & Dodds 1994). Not surprisingly 
they are also in accord with other measures of "shape" 
such as the zero-crossing of the cluster correlation function 
({(fee) = for Tec ~ 50/t-^Mpc, (Klypin & Rhee 1994)) and 
are reasonable fits to the APM angular correlation function 
w{9). For h < 0.5 we can reduce ag to acceptable levels, 
though as can be seen in Figure 1 we find it difficult to 
obtain ag < 0.7 without the inclusion of tensors. 

The amplitude of the mass fluctuations on l/t~^Mpc 
scales is consistent with estimates from the abundances 
of quasar and damped Lya systems at high redshift, 
very roughly fi > 3 today (Haehnelt 1993; Ma & 
Bertschinger 1994; Kauffmann & Chariot 1994; Mo & 
Miralda-Escude 1994; Subramanian & Padmanabhan 1994; 
Liddle & Lyth 1995). These are rather soft numbers to aim 
at, but the models we consider have enough small scale 
power to form objects early, in contrast with models like 
MDM where small scale power in the hot component is ex- 
ponentially damped. 

Figure 1 or Table 2 shows that a factor of 3-4 reduction 
in J2(l/t~^Mpc) over the standard CDM number is possi- 
ble, though it is difficult to obtain more suppression than 
this and still get high ai. (The values of J3(l /t~"'Mpc) show 
a similar dependence on the input parameters, so if J2 is in 
the acceptable range, we expect J3 will be also.) The "hot" 
velocity fields on Mpc scales are potentially the biggest prob- 
lem for these models, in common with most other Qo = 1 
models. In order to address the question of possible velocity 
bias (Carlberg & Couchman 1989), and to study this prob- 
lem in more detail, we are investigating the velocity fields 
for these models in cosmological N-body simulations. 



5 CONCLUSIONS 

We propose that the two models in Table 2 are currently the 
'best-buy' versions of CDM, and contend that these models, 
in which the parameters are varied by ~ 10% from their 
canonical values, should be the standard by which CDM is 
judged. Before abandoning the appealing simplicity of stan- 
dard CDM in favour of A or massive neutrinos (or more 
extreme modifications, such as isocurvature fluctuations or 
topological defects), it is these inflationary-inspired CDM 
models which should be confronted with observational data. 

Despite our belief that these mild variants of CDM are 
a reasonable fit to all currently available data, there are cer- 
tainly future tests which could definitively rule them out. 
The age determination is the one classical test that is rel- 
evant because it is largely the recent flurry of activity in 
Hubble constant determinations that is forcing reconsidera- 
tion of low-Q universes: if measurements of h become firm 
at values like 0.7-0.8, then there is no way to save CDM 
without at least introducing a cosmological constant. How- 
ever, conflicting claims abound, and we are presently reluc- 
tant to place any more emphasis on h = 0.8 ± 0.08 than 
A = 0.5 ± 0.05. 

The present situation on large scale flows is ambigu- 
ous (Dekel 1994; Strauss & Willick 1995; Shaya, Peebles 



& TuUy 1995), but improved data should yield a definitive 
measure of /? = flg'^ag. Given a reliable determination of ag, 
the degeneracy between Qo and ag could be broken. If indi- 
cations that Qo < 1 from these considerations are ultimately 
confirmed, then we will have to give up the aesthetic appeal 
of an Einstein-de Sitter universe. Another strong test will 
be the so-called 'baryon catastrophe' in clusters (White et 
al. 1993): if the fraction of baryons in clusters implies that 
Q_b/Qo — 0.2 then this is difficult for all high Qo models. 
Here one currently runs into uncertainty in the determina- 
tion of cluster masses as measured by the X-ray emitting 
gas. Lensing studies of clusters will be an important check 
on possible systematic errors in the available mass estimates. 

From the point of view of the microwave background, 
the most damaging observations would be: if the acoustic 
peak became a firm detection at smaller angular scales, since 
it is not possible to achieve this with a flat geometry; or if 
the peak has height ^ 4, which is hard to achieve with ra < 1 
and no cosmological constant. Another clear test which may 
perhaps be performed in the next few years is a more defini- 
tive measurement of the matter power spectrum shape. If 
the break in the power spectrum is really at larger scales 
than predicted by CDM, mild variations of the CDM hy- 
pothesis no longer suffice, and a more drastic revision of our 
current best-fit model for large-scale structure is required. 
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Coefficient of 

Quantity 1 h n hn 'h? h'h? 

Ao.05/<78 0.22 -0.28 0.17 -0.16 

(78 1-33 3.31 -2.40 -6.27 — — 

(71 11.1 31.6 -34.5 -79.9 41.4 — 

J2(l/i-lMpc) 205 1060 -738 -2899 536 — 

(78 1.33 3.19 - 5.54 -9.2 5 9.58 — 

(71 11.1 30.3 -55.0 -100. 103. — 

J2(l/i-lMpc) 203 1034 -1493 -5753 3355 8971 

Table 1. Fitting formulae for measures of power as a function of 
/i = /i — 1/2 and ri = 1 — ri, valid to 5% over the range 0.45 < 
h < 0.55 and 0.8 < ?i < 1. The error on the normalization from 
GOBE is 7.5% for (78 and ai and 15% for J2(l /i-^Mpc). The 
fits are divided into two cases: without (upper) and with (lower) 
tensors. The shape ratio is independent of the tensor contribution. 



Quantity 


"Obs" 


C,^=0 


C,^>0 


"CDM" 




327 ± 82 


282 


254 


356 


Ao.os/fs 


0.27± 0.02 


0.27 


0.25 


0.22 




0.5-0.8 


0.75 


0.77 


1.34 




> 3-4 


5.0 


5.6 


11 


^rms 


< 600 


740 


775 


1400 


-D200/-D10 


2-10 


4.0 


3.3 


5.5 



Table 2. Power spectrum measures for our preferred CDM mod- 
els: h = 0.45, n = 0.8 (without tensors) and h = 0.45, n = 0.9 
(with tensors). The values of Wrms (kms~-^) come from J2 and 
Eq. (1), assuming rieif = —1. The "Observational" numbers come 
from papers cited in the text. 
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1023 

White S.D.M., Navarro J.F., Evrard A.E., Frenk C.S., 1993, Na- 
ture, 366, 429 
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Zurek W.H., Quinn P.J., Salmon J.K., Warren M.S., 1994, ApJ, 
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Figure 2. The matter power spectrafor our CDM models (solid), 
along with an MDM model (short dashed) and ACDM model 
(long dashed) for comparison. The ACDM model has been multi- 
plied by to allow comparison with the data points (from Pea- 
cock & Dodds (1994)). The non-linear P(A;) for the CDM models 
is shown dotted. 



Figure 1. Contours of (Tg (solid), A(0.05)/(T8 (dot-dashed), g\ 
(dotted) and J2(l/i~^Mpc) (dashed) for models with (panel b) 
and without (panel a) tensors. The contour of (7\ has been sup- 
pressed in panel (b) for ease of viewing. The shaded region shows 
(Tg > 0.8 and the "standard" Cold Dark Matter model is marked 
by an asterisk. Our best fit model has (a) h ~ 0.45, n ~ 0.8, 
jC^ = and (b) h ~ 0.45, n ~ 0.9, /C^' = 0.7. 
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Figure 3. The radiation power spectrafor our CDM models with 
(solid) and without (dashed) a tensor component, normalized to 
the 2-year COBE data. 



